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I. INTRODUCTION 

The problem of finding the equations of motion of a two-body system subject to gravita- 
tional interaction has been intensively studied since the advent of general relativity. Because 
no exact solution is known, different frameworks and approximations have been developed 
out in the past: the present work lies within the post-Newtonian (PN) approximation of 
general relativity, see e.g. pQ and [2] for reviews and current situation of this approach. 

In particular the computation of the effective two-body dynamics within the PN approx- 
imation has been the subject of intensive research in the last three decades, both in the 
conservative and in the dissipative sector. At present the Arnowitt-Deser-Misner Hamilto- 
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nian ruling the conservative dynamics of a gravitationally bound, spin-less binary system, 
has been computed up to third PN order in [3], whose calculations have been later finalized 
and confirmed by [IH5] in harmonic coordinates, and also by [7J , where a different resolution 
of the source singularity has been adopted. 

From the 2.5PN order onwards, the dynamics of binary systems is modified by dissipative 
effects [Sj, associated with gravitational radiation. At present the energy flux emitted by a 
spin-less binary system has been computed up to 3.5PN order with respect to the leading 
order contribution [5]. 

The main result of this paper is the re-computation of the two-body effective action for 
the conservative dynamics at 3PN order for spin-less objects via an automatized algorithm, 
making use of a different method than those implemented so far to obtain such a result: 
we performed the computation within the framework of the Effective Field Theory methods 
(EFT) for non relativistic particles introduced in [10]. This method has already been applied 
with success on the conservative sector of the theory to determine the two-body 1PN [TU] 
and 2PN [TTJ Hamiltonian (see also (T2J for some investigations beyond 2PN), as well as to 
the study of the n-body 2PN [13] Hamiltonian (for 3-body systems, first obtained in [E]). 
On the radiative side, [15] showed an application of EFT methods to the calculation of some 
terms of the gravitational wave energy flux for spin-less systems up to 3PN order, and in [16] 
it is shown how present gravitational wave observations can set bounds on the fundamental 
parameters of the theory. The EFT approach has lead not only to re-derivation of old and 
established results, but also contributed to the computation of new spin effects in the 3PN 
conservative p2] (also obtained in |18j with more traditional methods) and non-conservative 
dynamics [19] . 

An EFT method possesses several features making it suitable for computing physical 
quantities in problems involving well-separated scales. This is the case for gravitationally 
bound two-body systems, where the different scales can be indentified in the size of the 
compact objects r s , the orbital separation r and the gravitational wave-length A with hier- 
archy r s < r ~ r s /v 2 < A ~ r/v, being v the relative velocity between the two bodies (we 
posit the speed of light c = 1). We neglect here radiation effects and consider point-like 
sources, comforted by the effacement principle [20] , which guarantees that finite size-effects 
come into play from 5PN order, for spin-less objects. The power-counting scheme outlined 
in JTU] enables to clearly separate scales and to frame divergences arising from zero-size 



3 



sources within a well under control field theory setting. Moreover EFT methods make sys- 
tematic use of Feynman diagrams and exhibits manifest power counting rules, enabling to 
implement PN computations within an algorithm automatizing the intricate computations 
involving hundreds of Feynman diagrams. Such an algorithm is completely general and it 
can be generalized to the calculation of higher-order dynamics. 

A key simplification enabled by EFT methods with respect to the standard ones is the use 
of perturbative techniques to integrate out the degrees of freedom responsible of mediating 
the conservative gravitational interaction (the so-called potential gravitons in the language 
of JTU]) to directly obtain the 2-body action, without the need to solve for the metric. 

This field of studies is nowadays of great phenomenological impact as the large interfer- 
ometers LIGO and Virgo have been operating until recently (fall 2010) at unprecedented 
sensitivity and a further increase in sensitivity is scheduled for their advanced runs, due 
to start in the year 2014, see e.g. [21] for a recent publication of coalescing binary signal 
search. As reported in [22J , reasonable astrophysical estimates make plausible the detection 
of gravitational wave signals from coalescing binaries by the advanced interferometers. 

In the search for signals from coalescing binaries the detector output is processed via 
standard match-filtering methods [23J, where the experimental data are confronted against 
banks of template waveforms. The result is particularly sensitive to the time varying phase 
of the oscillating gravitational wave signal [23]. In order to compute such a phase with 0(1) 
accuracy it is necessary to take into account PN corrections to both the energy and the 
emitted energy flux up to 3PN order at least. 

The frequency region accessible to the advanced versions of LIGO/Virgo (roughly from 
10Hz to 10 kHz) corresponds to the last stages (whose duration is at most few tens of 
minutes) of the coalescence of astrophysical objects like black holes and neutron stars and 
by then the orbit will have circularized, as eccentricity decay faster than orbital separation 
[21)] . enabling the the energy and flux functions to depend on the single parameter v. 

The inspiral phase of the coalescence is bound to end when the two constituents of the 
binary system merge into a single object and the perturbative PN analysis cannot be followed 
beyond the merger. Here however numerical techniques have been proven useful, see e.g. 
[26] for reviews, to construct complete analytical waveforms from inspiral to merger and the 
subsequent ring- down, as it is conventionally named the last phase in which the final object 
produced by the merger undergoes damped oscillation before settling down to quietness. 
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The ring-down phase also admits an analytical perturbative description 

The present work represents a decisive step towards the computation of the Hamiltonian 
dynamics at fourth PN order, which is a key ingredient in the construction of template 
waveforms to be used in actual experimental searches, like the the Effective One Body 
waveforms |28] which have been used in the recent 1291. 



The paper is organized as follows. In sec. II A we give an overview of the effective field 



theory methods for gravity, which are applied in sec. Ill to the determination of the 2-body 



effective action at order 3PN. We summarize and conclude in sec. IIVI 

II. MAIN 
A. Lagrangian and Feynman rules 

The relevant physical scales for the 2-body problem in gravity are the total mass M, 
the separation r and the relative velocity v. By nPN correction it is conventionally meant 
corrections of order v 2n ~ (G^M /r) n , where Gn is the standard Newton constant and the 
virial relation has been applied, showing that an expansion in v 2 is at the same time an 
expansion in the strength of the gravitational field. 

At 3PN order, the EFT description of massive compact objects in binary systems takes 
them as non-dynamical, background point-like sources: quantitatively this corresponds to 
particle worldlines interacting with gravitons. The action S we consider is then given by 

S = Seh + Sqf + S pp , (1) 

where the first and the third terms are, respectively, the <i-dimensional Einstein-Hilbert 
action and the worldline point particle action 

Seh = 2A 2 y d d+1 x^9 R(g) , 

i=l,2 J i=l,2 J 

with A~ 2 = 32ttG, being G the (^-dimensional gravitational constant [37] : as dimensional 
regularization will be needed, it is necessary to keep d generic until the actual value of 
Feynman integrals will have been computed. The goal is to compute the effective action S e ff 
for particles alone, with gravitons mediating interactions being integrated out by standard 
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perturbative methods, i.e. by the aid of Feynman diagrams. For the gauge fixing term, we 
follow PI]: 



S 



GF 



-A* 



X\ 



-<7r„r" 



(3) 



^a/39 a i corresponding to the same harmonic gauge adopted in pQ. Still following 



with 

[llj . we adopt the standard Kaluza-Klein (KK) parametrization of the metric [3U] (a some- 
how similar parametrization was first applied within the framework of a PN calculation in 
[3~T]). suitable for a nonrelativistic expansion around a Minkowski metric: 



with = 6ij + o-y/A, c d 
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and z,j running over the d spatial dimensions. 



In terms of the metric parametrization Q, each world-line coupling to the gravitational 
degrees of freedom (ft, reads 
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(5) 



and its Taylor expansion provides the various particle-gravity vertices of the EFT. 

Also the pure gravity sector Sbuik = Seh + Sgf can be explicitly written in terms of the 
KK variables, Seh m KK variable has been first derived in [3D]. We report here only a part 
of Sbuik, including all the terms that are needed for the 3PN calculation: 
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e a 



(6) 



where = Aj } i — Aij and indices are raised and contracted by means the <i-dimensional 
metric tensor 7 (although in most terms one can just use 5^ because the neglected parts 
are not needed at 3PN)[38]. All spatial derivatives are meant to be simple (not covariant) 
[39] and, when ambiguities might raise, gradients are always meant to act on contravariant 
fields (so that, for instance, VA = 7 y A,j and F^ = y k 'y jl F i jF k i). 

All the Feynman rules needed for our calculation can be read off eq. ([6]). In particular: 
the first line gives the a and (f) propagators, as well as all the <p n a m interaction vertices; the 
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second line gives the A propagator and, together with the third, contains enough information 
to compute all the other vertices needed for the 3PN calculation. 

Borrowing the definitions from [10], the two bodies exchange potential gravitons, respon- 
sible for binding the system as they mediate instantaneous interactions: their characteristic 
four-momentum k^ scales thus as ~ (k° ~ v/r, k ~ r _1 ) and we neglect altogether the 
emission of gravitational waves because we are interested in the conservative part of the 
binary system dynamics. When a compact object emits a single graviton, momentum is ac- 
tually not conserved and the non-relativistic particle recoils of a fractional amount roughly 
given by — |&|/H — h/J, where J ~ mvr is the angular momentum of the system: 

however, for macroscopic systems such a quantity is negligibly small. 

In order to allow manifest scaling it is necessary to work with the space-Fourier trans- 
formed fields 



k (t) = J d d x0(t,x)e- ikx , 
A i>k (t) = J d d ^,(i,x)e- ikx , (7) 

a ij>k (t) ee J d d iff 1J (t,x)e- !kx . 

The fields defined above are the fundamental variables in terms of which the Feynman 
graphs are going to be constructed; the action governing their dynamics can be found from 
eqns. (5|6) by direct substitution and replacement of spatial derivatives by the appropriate 



ik factors. 

From the quadratic part of the Lagrangian the propagators can be written as: 

1 



P[Mta)Mtb)} = 
P[Ai tk (t a )A jtV (t b )] = ^Sij 

PWij,'k{ta) 'kl,'k?(tb)] = ~^PijM 



} x (2n) d 6 d {k + k')P(k 2 , t a , t b )5(t a - t b ) , (8) 



2 

where P ijM = - [8^8^ + 5u5 jk + (2 - c d )<%<!> fc /] and 

is the full relativistic propagator, which can be thought as an instantaneous nonrelativistic 
part plus insertion terms involving time derivatives. Only few of these time derivative terms 
need to be included in a Feynman diagram at a given PN order. 
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B. Feynman diagrams 

We have exposed all the ingredients needed to obtain the 2-body effective action S e ff 
with manifest power counting in G and v at the third post-Newtonian order. This can be 
done by integrating out the graviton fields from the full action derived above 

S e ff = J D(pDaijDA k exp[i(S EH + S GF + S pp )] . (10) 

We can now proceed to lay down all the Feynman diagrams relevant for our computation, 
then by applying the Feynman rules derived from the above interactions the amplitudes will 
be computed, collecting all results belonging to the same PN order. 

1. Topologies 

Following the procedure and the terminology of [11] , we first select the topologies having 
the appropriate power of G for being potentially relevant from a given PN order on, that is 
Qi n + l ) f or rith-PN order. The attribution of a power of G to a topology is made according 
to the following rule: each vertex involving n gravitational fields gives a factor G n/2 ~ l if it 
is a bulk one, and a factor G n l 2 if it is attached to an external particle. 

Concretely, there is just one topology contributing at the Newtonian order G: the one 
representing the single graviton exchange. At 0(G 2 ) order two other topologies (modulo 
switch of the external particles) enter the game but actually just one of these contributes 
at 1PN since the other involves an extra v 2 factor, which is equivalent to increasing by one 
unity the PN order. As underlined in [TTJ, this is a very welcome effect of using the KK 
parametrization, also with respect to the Arnowitt-Deser-Misner one, see |32j . 

At 2PN order, 5 new G 3 topologies have to be added, while 4 additional ones come into 
play only from 3PN on because of v 2 suppressions. All these 12 topologies contributing to 
G 3 order have been already displayed in [11] , and they are also reported here in figure [TJ 

To complete the 3PN order the introduction of G A topologies is required. The generation 
of the new topologies at any given order n is done iteratively from the ones at the previous 
order n — 1: to any G n_1 topology a new propagator is attached in all possible ways, provided 
that one extremum ends on an external particle, and the second extremum is placed in one 
of the following three locations: i) an internal vertex, ii) a vertex attached to the other 
external particle, iii) a propagator, so to "break" it into two propagators ending into the 
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Figure 1: The 12 topologies up to G 3 as they are produced by a Mathematica code. The last 4 
start contributing at 3PN order. Here and in the following figures, we count them form left to 
right, from top to bottom. 

same newly-born 3-legged internal vertex; barring in all cases the possibility that internal 
graviton loops are created, as they would give rise to uninteresting quantum corrections (i.e. 
proportional to h/ J <C 1). 

In this way one actually obtains a redundant set of G n topologies because many of them 
are equivalent, i.e. they are skeletons of the same Feynman graphs, then the algorithm to 
eliminate equivalent graphs can be straightforwardly implemented on a computer program. 

We have indeed implemented the whole procedure into a Mathematica code [33] and 32 
G A topologies have been found, only 8 of which turned out to be relevant at the 3PN level 
of accuracy, as the remaining 24 contribute at higher PN-orders because of v 2 suppressions. 
Moreover, as figure [2] shows, none of the 8 3PN-relevant topologies is really new, since they 
are actually factorisable in lower order sub-topologies. We have also worked out, for the 
sake of future developments, the situation at 4PN and summarized the number of topologies 
in the left part of table [1} 

That makes a total of 20 topologies that must be taken into account for the 3PN calcu- 
lation, and only 2 of them (the two G 3 ones involving 4-legged vertices, the last two ones of 
fig. [T]) are genuinely new ones, i.e. non-factorisable into subtopologies. 

2. Graphs 

The next step is to generate Feynman graphs from all topologies by replacing all the 
not-yet-specified propagators by the proper 0, A or a ones. This is the point where powers 




Figure 2: The 32 G 4 topologies. Only the first 8 contribute at 3PN order, and they are all 
factorisable into simpler subtopologies. 
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V V 
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1 
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G 2 




1 1 


G 2 


1 


8 


7 7 


G 3 




5 4 


G 3 


5 


48 


159 .. . 


G 4 




8 21 


G 4 


8 


299 




G 5 




50 


G 4 


50 







Table I: On the left: number of topologies entering at a given PN order. On the right: number 
of diagrams that start to contribute to the effective action at a given G and v power. The 63 
diagrams whose contribution starts at 3PN are written in boldface characters. 

of v (traced by time derivatives in the Lagrangian) enter the game, both in external and bulk 
vertices and in the relativistic propagators. The lowest power of v 2 appearing in a graph, 
together with its G-order, determines the lowest PN order at which such a graph must be 
taken into account. As different vertices carry different time derivative structures, see for 
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instance eq. (J6J), topologies of the same G-order, or even the same topology, can generate 
Feynman graphs whose lowest PN orders differ. In particular some topologies allow only 
graphs containing v 2 corrections at least, thus entering the game at higher PN order than one 
would have inferred by a naive G-power counting; this is the origin of the above mentioned 
simplification introduced by the use of KK variables, which reduces the number of topologies 
(and thus of graphs) to be taken into account at any given PN order. This welcome feature 
becomes more and more important at higher PN orders. 

The generation and classification of all the Feynman graphs needed at a given PN level is 
an automatizable procedure, as it suffices to try all possibilities (although one can envisage 
shortcuts towards more efficient algorithms) and find out the lowest power of v 2 of each 
graph using the Feynman rules implicit in the Lagrangians 
has thus generated all 80 the diagrams giving a contribution at 3PN (as well as other 515 
for the 4PN case): 63 of them are new ones, while the other 17 have been computed at 2PN 
order by Gilmore and Ross [TlJ[40j. The right part of table [T] shows the classification of the 
new diagrams according to their G and lowest v 2 powers. 



5p). Our Mathematica code 



III. THE RESULTS 

We use a Mathematica [33] code together with EinS [34J and Feyncalc |35J softwares to 
generate and to compute all the diagrams. Here we show the results concerning the 3PN 
part only |41j . 



A. 0{G) diagrams 

They are shown in figure [3] and they can all be done by means of the dimensionally 
regularized Fourier formula 



ipr 1 T(d/2-a) /2 X d ~' 2n 



P P 2a (An) d/2 I» 



(11) 



(with j p = J 7§J^j and r = xi — x 2 ), as well as its generalization to the case where factors 
of p M 's appear in the integrand; the corresponding formulae can be obtained by taking the 



appropriate number of r-gradients of equation (11) 



The result is then multiplied by a symmetry factor accounting for: 
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2 




Figure 3: The three 0(G) diagrams. The (f), A and a propagators are represented respectively by 
blue dashed, red dotted and green solid lines. From the left to the right, they contain f 6 , i> 4 and 
v 2 corrections under the form of vertices expansion and/or double time derivative insertions in the 
propagators. The number in the right part of each diagram is its multiplicity factor. 



factorials coming from the development of the exponential function contained in (10); 



• a factor 1/2 for diagrams symmetric under particles exchange, in order not to over- 
count when performing symmetrization (i.e. + 1 -H- 2) in the very final stage of the 
calculation; 

• the number of possible Wick contractions giving rise to the same diagram, with the 
convention that equivalent legs attached to an internal vertex are indistinguishable 
(because we chose to incorporate the corresponding factors already in the internal 
vertices definitions). 

The results of the integration have been evaluated at d — 3 and multiplied by the imaginary 
unit i to give the following potential terms for the 3PN effective Lagrangian (G^ being the 
value of the fe-th G a diagram, listed as in figures [MoT) : 
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G N m x m 2 f a x • a 2 + 3n.axn.a 2 4 



G 1 = — — | — : — r + [ai.a 2 r.a 2 - 2ai.a 2 r.ai - a 2 .air.vi - a 2 .air.v 2 

— 2r.aiVi.a 2 — 18vi.aiVi.a 2 — 2ai.a 2 t> 2 + %a\v\ — 2ai.a 2 Vi.v 2 — vi.a 2 v 2 .ai — 14r.aiv 2 .a 2 
-lllvi.aiv 2 .a 2 - %a\v\ - 7ai.a 2 i» 2 ] r 2 + 11v\ + 2(vi.v 2 ) 3 - 2ai.a 2 (r.vi) 2 - ai.a 2 (r.v 2 ) 2 
+16w 2 (vi.v 2 ) 2 — (r.ai) 2 r.a 2 — r.a 2 r.air.vi — r.a 2 r.air.v 2 — 2ai.a 2 r.vir.v 2 
— 7r.air.viVi.a 2 — 3r.air.v 2 vi.a 2 — 3(r.a 2 ) 2 t> \ — 2r.air.a 2 t> 2 + 88r.v 2 vi.aiv 2 
+22f^V!.v 2 — 2r.a 1 r.a 2 v 1 .v 2 + 32r.v 2 v 1 .a 1 v 1 .v 2 + 16r.v 2 f 2 v 2 .a x + 8r.v 2 v 1 .v 2 v 2 .a 1 
— 41r.air.viv 2 .a 2 — 17r.air.v 2 v 2 .a 2 + 3(r.ai) v 2 + 3Av 1 v 2 — 7r.air.a 2 t> 2 + 124r.v 2 vi.ait> 2 

+23f 2 vi.v 2 f 2 + 20r.v 2 v 2 .aiW2 + 2r.air.a 2 (n.vi) 2 + r.air.a 2 (n.v 2 ) 2 + 8(n.v 2 ) 2 t> 4 
— 22n.vin.v 2 t> 4 — 6n.vin.v 2 (vi.v 2 ) 2 + 2r.air.a 2 n.vin.v 2 — 4r.ai(n.v 2 ) 2 Vi.v 2 
+12(n.v 2 ) 2 f 2 Vx.v 2 — 16n.v 1 n.v 2 v 2 v 1 .v 2 — 8(n.v 2 ) 3 v 2 .a!r/3 — 4r.v 1 (n.v 2 ) 2 v 2 .a 1 
— 20r.ai(n.vi) 2 f 2 — 10r.ai(n.v 2 ) 2 v 2 — 16r.ain.vin.v 2 V2 — 20(n.vi) 2 v 2 V2 
— 39n.vin.v 2 f \v 2 + 2r.ai(n.v 2 ) 4 + 4r.ain.vi(n.v 2 ) 3 + 20(n.v 2 ) 4 t> 2 + 4n.vi(n.v 2 ) 3 t> \ 

+9(n.v 1 ) 2 (n.v 2 ) 2 v 1 .v 2 - 5(n. Vl ) 3 (n.v 2 ) 3 } , 

! G N m 1 m 2 ( 4 ai.a 2 . 

G 2 = < r — h [ai.a 2 r.a 2 - 2a 2 .a!r.vi - 2vx.a 2 r.ai - 2a 2 .aii\v 2 

4r [ 3 

— 3ai.a 2 vi.v 2 — 2v 2 .aiVi.a 2 — 14vi.aiVi.a 2 — 3t> 2 ai.a 2 — 5f 2 ai.a 2 — 3ai.a 2 r.ai] r 2 
— 2a 1 .a 2 r.v 1 r.v 2 — 2r.a 1 v 1 .a 2 r.v 2 — r.a!r.a 2 V!.v 2 — 3a 1 .a 2 (r.v 1 ) 2 — 6r.a 1 v 1 .a 2 r.v 1 
— ai.a 2 (r.v 2 ) 2 + 12vi.air.v 2 vi.v 2 + 6w 2 v 2 .air.v 2 + 8v 2 .air.v 2 vi.v 2 + 10w 2 v 2 .air.v 2 
— 2r.ai(n.v 2 ) 2 vi.v 2 — 2v 2 .air.vi(n.v 2 ) 2 — 2rv 2 .ai(n.v 2 ) 3 + 2(vi.v 2 ) 3 + 6w 2 (vi.v 2 ) 2 

+6w 4 Vi.v 2 + w 2 -u 2 V!.v 2 - 4n.V!n.v 2 (vi.v 2 ) 2 - 6w 2 n.V!n.v 2 Vi.v 2 + 3(n.Vi) 2 (n.v 2 ) 2 Vi.v 2 } , 

i G^mim 2 f 2 
G 3 = |[2ai.a 2 vi.v 2 + 2v 2 .aiVi.a 2 — 4vi.aiv 2 .a 2 J r 

-4v 2 .air.v 2 Vi.v 2 + 4v 2 Vi.air.v 2 - (vi.v 2 ) 3 - 2-y 2 (vi.v 2 ) 2 + 2v^vj + v 2 vlv 1 .v 2 

+n.vin.v 2 (vi.v 2 ) 2 — vjvln.v 1 n.v 2 } , 

where n = r/r and the Newton's G^ constant is related to the d— dimensional gravitational 
constant G by the following relation involving the arbitrary subtraction lenght scale L: 

G = G N L d ~ 3 , (12) 

which gives trivially G — > Gn in this case, but gives also a nontrivial contribution in presence 
of divergences, see next section. 



13 




Figure 4: The 16 0(G 2 ) diagrams. Diagrams from 2 to 9 actually start contributing at 2PN 
order, and must be consequently evaluated at their next-to leading order in v 2 (which involves an 
appropriate vertices expansion or one double time derivative insertion in the propagators) to catch 
their 3PN contribution. Analogously, the first diagram enters already at 1PN order, so it must be 
expanded up to u 4 . As in fig. [3J the (f), A and o propagators are represented respectively by blue 
dashed, red dotted and green solid lines and the number in the right part of each diagram is its 
multiplicity factor. 

B. 0(G 2 ) diagrams 

The first, second and tenth diagram of figure [4] can be factorized in terms of lowest order 
ones, all the others need the use of the well known formula 

1 r(d/2 - P)T{d/2 - q)r(q +J3- d/2) ^ 2a _ 2g ^ 



/ Pl V\ a iv ~ Pi) 2/3 (4nf 2 r(a)rG0)r(d -a-fi) 

and of its generalization to inclusion of p^'s factors in the integrand (up to four different 
indices are needed at 3PN), which can be worked out by means of appropriate contractions 
with p M 's. The result is: 

G\ = "J" 1 ™ 2 [2r 2 ai.a 2 + 36v 2 r. ai + 8 Vl .a 2 r. Vl - 8r.a 2 (n. Vl ) 2 - 8( Vl .v 2 ) 2 - 24v 2 Vl .v 2 

+lUvlvl - 32v* + 49^2 + 48v 2 n.v 1 n.v 2 - 24(n.Vi) 2 Vi.v 2 + 32n.Vin.v 2 Vi.v 2 



-28^(n.vi) 2 - 24w 2 (n.vi) 2 - 32(n.vi) 2 (n.v 2 ) 2 + 32(n.vi) 3 n.v 2 - 8(n.vi) 4 ] , 
G 2 = [r ai.a 2 + 2vi.a 2 r.vi — 2(vi.v 2 ) — 4f 1 v i .v 2 + 3u 2 Vi.v 2 — 2(n.vi) vi.v 2 



-4n.Vin.v 2 Vi.v 2 ] 
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24r 2 log (j-^J (2ai.a 2 + aj) - r 2 (32ai.a 2 + 22a 



G 2 N m\m 2 lnn 2 G^mfo 

~3(d-3) [2ai - a2 + ai J + ^^ 
+48r.air.a 2 — 180t> 2 r.ai — 96v 2 .air.v 2 — 156i> 2 r.ai — 12(r.ai) 2 + 96r.ai(n.v 2 ) 2 

+54vi.a 2 r.vi + 288w 2 r.a 2 - 54r.a 2 (n.vi) 2 - 54(vi.v 2 ) 2 + 150v 2 Vi.v 2 - 369^1 - 78v- 

+576f 1 (n.v 2 ) — 300f 1 n.vin.v 2 + 42(n.Vi) Vi.v 2 + 216n.vin.v 2 vi.v 2 + 243t> 2 (n.vi) 

+222w 2 (n. Vl ) 2 - 216(n. Vl ) 2 (n.v 2 ) 2 - 56(n. Vl ) 3 n.v 2 + 20(n. Vl ) 4 ] , 
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G\ = — ^ I — - Tr 2 ai.a 2 + 3r.air.a 2 + 12r.aiVi.v 2 — 3v 2 .air.v 2 — 3wor.ai + 6r.ai(n.v 2 ) 2 

br z 

+llvi.a 2 r.vi + 25-u 2 r.a 2 - 2r.a 2 (n.vi) 2 - ll(vi.v 2 ) 2 + 6v 2 vi.v 2 - 25vfvl 
— 3t>2Vi.v 2 + 50f 2 (n.v 2 ) 2 — 24f 2 n.vin.v 2 — 12(n.Vi) 2 Vi.v 2 + 26n.Vin.v 2 Vi.v 2 
+2t>2(n. vi) 2 + 3v 2i1.V1n.V2 — 8(n.vi) 2 (n.v 2 ) 2 — 4(n.vi) 3 n.v 2 ] , 
Gl = Gn ™^ 2 [lOv 2 .air.v 2 + 2w 2 r.ai - 4r.ai(n.v 2 ) 2 + 6(vi.v 2 ) 2 + l%v\v 2 2 + 3v 4 
-24-y 2 (n.v 2 ) 2 - 12n.vin.v 2 vi.v 2 - 3-u 2 (n.v 2 ) 2 + 12(n.Vi) 2 (n.v 2 ) 2 ] , 

G l = - GN ^f 1/2 V 2 ( a i- a 2 + 2a 2 ) + 27w 2 r. ai - 9^r.ai + 14 Vl .a 2 r. Vl + 4w 2 r.a 2 

-14r.a 2 (n.V!) 2 - 14(vi.v 2 ) 2 - 39t> 2 V!.v 2 - 22v\vl + 9vf + 8w 2 (n.v 2 ) 2 + 78v 2 n.Vin.v 2 

+56n.vin.v 2 vi.v 2 + SO^fn.Vi) 2 - 24w 2 (n.vi) 2 - 56(n.vi) 2 (n.v 2 ) 2 + 8(n.vi) 4 ] , 

^ 2 8G 2 N m 2 m 2 , AG 2 N m\m 2 
G 7 = 1 a a i a 2 + - 2 



..2. ' - • r 



d - 3 — - ' f ai ' a2 I' " 41 ° g UJ ) " 2r ' aiVl ' V2 

+4v 2 .air.v 2 — 4vi.a 2 r.vi — v 2 r.a 2 + 4(vi.v 2 ) 2 + 3v 2 vi.v 2 + v\v\ + t>2Vi.v 2 

— 2v 1 (n.v 2 ) + 2t> 1 n.vin.v 2 + 2(n.vi) vi.v 2 — 8n.Vin.v 2 Vi.v 2 J , 

Gl = Gn ™ 1 ™ 2 Tg^r.a! - 2k; 2 vi.v 2 - 27v 2 v% - 33v 4 + 42w 2 n.vin.v 2 + 96(n.vi) 2 vi.v 2 

9r 2 

+54w 2 (n. Vl ) 2 + 24w 1 2 (n.v 1 ) 2 - 128(n. Vl ) 3 n.v 2 + 56(n. Vl ) 4 ] , 



i2 



AG 2 N m 2 l m 2 2 2G 2 N m\m 2 

69 " d-3 ai + 7 2 



4r 2 a 2 ( 1 - log ( ^- ) ) + 8w 2 r.ai - 6v 2 Vi.v 2 



-3v\v 2 2 + v\ + 12w 2 n.vin.v 2 - 9w 2 (n.vi) 2 ] , 

2 QG 2 N m\m 2 r , , 2 2 21 
G10 = ^ L(vi.v 3 ) - v lV2 \ , 

2 SG^m^ma 
GJi = 7, v 2 .a!r.v 2 , 



r 



G12 = — - [vi.a 2 r.vi - (vi.v 2 ) 2 + -y 2 Vi.v 2 - 2(n.Vi) 2 Vi.v 2 + 2n.Vin.v 2 Vi.v 2 ] , 

G13 = — ^~2~^ — [4r.aiVi.v 2 + 6vi.a 2 r.vi — f 2 r.a 2 — 6(vi.v 2 ) 2 + 2v 2 vi.v 2 + t> \v\ 
— 2v 2 (n.v 2 ) 2 + f 2 n.vin.v 2 — 4(n.Vi) 2 Vi.v 2 + 12n.Vin.v 2 Vi.v 2 ] , 

Gl = 2 ^>1^ [v 2 v 2 _ 2(vi Va) 2 + V 2 (n . va )2] , 

G 2 5 = - 8 ^^, 2 r.a 1 , 

G 2 6 = ^^K-2(n. Vl ) 2 ]v,v 2 . 

Divergences appear at this order from the terms G|, G 2 and G 2 ,, under the shape of (d — 3) 
simple poles, giving rise also to logarithmic factors involving the subtraction scale L = 



16 



Le 7 / 2 /v^4vr , where 7 = .577216 ... is the Euler-Mascheroni constant and L is the length 



scale written in eq. (12). Following [Ij, in sec. HIE we will get rid of the poles by means of 
a shift of the worldline parameters; as to the logarithms involving the arbitrary subtraction 
scale, they are known to remain in the Lagrangian and also in the equation of motion for r, 
which is however a gauge-dependent variable. They however drop out in the expression of 
observables, like the energy E of the system as a function of the orbital frequency u, which 
are gauge-independent quantities as both E and u can be measured asymptotically far from 
the binary system. 



C. C(G 3 ) diagrams 



Most of the diagrams of figure [5] can be computed via recursive application of equation 



(13), with the exception of the 5th diagram and all the ones sharing the same H-shaped 



topology. In this case one has to work on the integrand by parts to obtain the following 
recursive relationship 

( 1)P2 Pi a (p - pi) 2 ^ 7 (p - p 2 ) 25 (pi - p 2 ) 2e 



= [ 1 (I(a-l,P, 1 + l,S,e)-I(a,P, 1 + l,5,e-l))+ ( 14 ) 
5(I(a,P-l, 1 ,6+l,e)-I(a,P, 1 ,8 + l,e-l)))/ (26 + 7 + 6 -d) , 

and its generalization to the inclusion of p± 2 factors in the integrand. The above formula 



(14) also holds exchanging in the right hand side <H- {7,5}) or ({a, 7} -H- {/3,5}) 



and it can be iterated to obtain an integral of the type eq.(13). 

As in the G 2 case, here also poles and logarithmic terms are obtained, which again are 



non physical, see sec. |III E| The output of our code for these diagrams is: 
Gl = ^^f 2 [l2v 2 - 3 Vl .v a + 7n. Vl n.v 2 - 4(n. Vl ) 2 ] , 
Gl = _ fiv m i m 2 V2 + 3y 2 + Qv 2 _ 3n Vin V2 + 2(n. Vl ) 2 ] , 



Gl = 3^j^p [lOvi.Va - 13vf - 30n. Vl n.v 2 + 39(n. Vl ) 2 ] 



+ 



G\ 



" N mim 2 



900r 3 



90 



log (j^J (I3v 2 - lOvi.Va + 30n.vin.va - 39(n.vi) 2 ) 



+ 850vi.v 2 - 2401w 2 - 450w 2 - 3450n. Vl n.v 2 + 3243(n. Vl ) 2 ] 

2 G%m\r" 2 
128r 3 



G\ = - ^ff* [ Vl , V2 + 2v 2 - 3n. Vl n.v 2 - 6(n. Vl ) 2 ] , 
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Figure 5: The 53 C(G 3 ) diagrams. The first 5 enter at 2PN order, and must be consequently 
evaluated at their next-to leading order in v 2 , which involve an appropriate vertices expansion 
or one double time derivative insertion in the propagators. As in fig. [3] and [4j the </>, A and a 
propagators are represented respectively by blue dashed, red dotted and green solid lines and the 
number in the right part of each diagram is its multiplicity factor. 

G 3 = ^^[ 9 ^ W _ ViV2 + 3nVinV2 _3 (nVi)2) 

+ 50vi.v 2 - 644v 2 - 294n. Vl n.v 2 + 348(n. Vl ) 2 ] , 

n% — ^7 _ ^li _ G\ 2 _ ^G\m\m\ 

^3 = _^23 _ 3^29 _ "iG\ % _ 2G%m\m 2 ^ 
8 /i 1« 32 r 3 1 2 ' 



Gj 




2 


4 


^23 


3^29 


4 


16 


^34 


G%m\i 


8 




3^28 _ 


3^45 


4 


8 



Gl = ^ = ^ _ (n . Vl) 2] 

G io " 4 " - g " r3 v i > 

2G 3 ^2_2 



G 3 3 = ^"™i m * [ Vl . Va _ 2n. Vl n.v 2 + (n. Vl ) 2 ] 
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8Gim?mn 



3 uk* N uu l uu 2 n< 2 

^14 - ^3 U l ' 



6G 3 2_2 



G\ b = -3GL = -^F^K-2(n. Vl ) 2 ] , 
G? 7 = -^I^H[v 1 .v 2 + ^-2n.v 1 n.v 2 -2(n.v 1 ) 2 ] , 



Gl s = 2 ^g^ [4(n.v 1 ) 2 -v 1 .v 2 -, 2 ] , 

y^(3 /^3 /^<3 /nr3 n 

^19 ~~ Lj 20 ~~ ^33 ~~ ^39 ~~ u > 

Gl = _ G W m i [ w 2 _ 4vi V2 + 14n . Vin . V2 _ 7(n.V!) 2 ] , 

GL = [v 2 _ 4vi>V2 + 12nvinV2 _ 3(n-Vi) 2] 

GL= ^^[r;?-2(n.v 1 ) 2 ], 
gL = 3G|^n, ^ _ (n va)2] 

o 2G%m%m2 r 2 , x21 

G 26 = -3 [vi.v 2 + «! - 3n.vin.v 2 - (n.vi) J , 

3 G 3 N mlm 2 

G 27 = -3 [V1.V2-n.V1Il.V2J , 

= K - 4vi.v 2 + 12n.vm.V2 - 5(n. Vl ) 2 ] , 

C31 = - ^3^3 ^ ~ 5Vl V2 + 15nv * nv 3 " 3 ( n - v i) 2 ] 

_ 2^m|m 2 3Q bg ^ (5vi ^ _ ^ _ 15n va + 3(n Vi)2) 

-25vi.v 2 - (vf) + 225n.vin.v 2 - 57(n.vi) 2 ] , 

^ 3 8G 3 N mjml r . 
Gg 5 = ^ [vi.v 2 - n.V!n.v 2 J , 

Gl 6 = [3tt 2 (vl - V1.V2 + 3n.vm.V2 - 3(n. Vl ) 2 ) 

+60vi.v 2 - QOvf - 148n.Vin.V2 + 148(n.Vi) 2 ] , 

G^ Tfl^Tfl^ 

^37 = ^77~3 — ~ ( Vl - V2 ~~ 3n.vin.v 2 ) - 32vi.v 2 + 32n.vin.v 2 ] , 



n 2 G 2 ! 2^2 



2r 3 

G 38 = " "j^™ 2 H - 3(n.vi) 2 ] , 
G 3 40 = ^ = ^.3 [vi.v 2 - 3n.vin.v 2 ] , 

G 3 2 = ^l^^^ + s^^^^^.^^^ 

G^* TYh^ TYb^ 

G 43 = 2 K ( v i- v 2 + v\ - 3n. Vl n.v 2 - 3(n. Vl ) 2 ) + 4vj.v 2 - Av\ + 4n. Vl n.v 2 - 4(n.vi) 2 ] , 
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G\ 4 



AG 3 N mlm 2 

(d-3)r 3 
2G 3 _3 



\y\ — vi.v 2 + 3n.Vin.V2 — 3(n.vi) 2 ] 



" N m\m 2 



6 log (jj^j ( v i- v 2 — vf- 3n.vin.v 2 + 3(n.V!) 2 ) 



r 3 



r 3 



^49 



+5u 2 — 5v!.v 2 + 21n.vxii.V2 — 21(n.V!) 2 ] , 

AG%m\m 2 
~3(d-3)r 3 

2G 3 v mfm 2 
+ 3^3 

-Vi.v 2 - v\ + 9n.vin.v 2 + 9(n.Vx) 2 ] 

2G%m 3 m 2 2G%m\m 2 
5(d-3)r 3 L 1 V ; J 
AG%m\m 2 



[vi.v 2 + v\ — 3n.vin.v 2 — 3(n.vi) 2 ] 
6 log (j^j ( v i + v i- v 2 — 3n.vin.v 2 — 3(n.vi) 2 ) 



25r 3 



17v'f - 56(n. Vl ) 2 - 15 log ( — ) (v\ - 3(n. Vl ) 2 ) 



15(d- 3)r 3 
r 



4G 3 v m 3 m 2 

[vi.v 2 — 3n.Vin.v 2 J [12vi.v 2 — 66n.vin.v 2 



75r 3 



15 log ( — ) (vi.v 2 — 3n.vin.v 2 ) 



r 3 



^51 



2G' 3 N m\m 2 
5(d- 3)r 3 
r 

2G' i N m\m 2 
15(d-3)r ; 



•15 log 



v x .v 2 — 3n.V!n.v 2 ] - 
V1.V2 — 3n.Vin.V2 



2G%m\m 2 
25T 3 



[17vx.v 2 — 56n.vxn.V2 



[vl - 3(n. Vl ) 2 ] + 2 °^ m2 [23vl + 36(n. Vl ) 



+ 15 log ( — ) (vf - 3(n. Vl ) 2 ) 



^52 



2G 3 N m\m 2 
3(d-3)r 3 
r 



[vi-v 2 - 3n.vin.v 2 ] 



2G 3 N m\m 2 
3T 3 



v x .v 2 + 6n.Vxn.v 2 



C53 



+3 log (vi.v 2 - 3n.vin.v 2 ) 

G ^ m?m2 [«? - 3(n. Vl ) 2 ] + °\™}™ 2 [I7vl - lOOvx.va - 81(n. Vl ) 2 + 300n. Vl n.v 2 



10(d- 3)r 3 300r 3 
- 901 °g(£;) K -3(n. Vl ) 2 ) 



D. C(G 4 ) diagrams 



The G 4 topology set shown in figure [6] does not present any new difficulty, and the 
integrations can be done straightforwardly: 
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Figure 6: The 8 0(G 4 ) diagrams. As in fig. [3j |4]and[5j the (f>, A and a propagators are represented 
respectively by blue dashed, red dotted and green solid lines. The number in the right part of each 
diagram is its multiplicity factor. 



G\ = G 



3G\ _ Gj _ 3GI _ G%m\m\ 
1 ~ 2 8 " 1 " ■><■' 



4 G\ _ G%mfm 2 



r .3 



Gi = 



24r 4 



E. Reproducing the 3PN Lagragian 

Adding up all the terms shown so far a Lagrangian containing poles and logarithmic terms 
is obtained. The poles can be removed by performing the shift of the worldline parameters as 
given by eq.(1.13) of [6J, acting over the source positions xi 2 according to xi — >■ x'-l = xi+£i 
with 



6 - yG 2 ™? 



1 , /r\ 327 

-2 log 



ai , (15) 



d-3 \L Q J 1540 

and analogously for £ 2 [42J. 

Finally, in order to get the exact form of the 3PN Lagrangian as it is written in eq. (174) 
of [1] (and originally derived in [5]), one has to act on factors like ai • a 2 , ai • a 2 and similar 



21 



by means of the standard double-zero technique, see e.g. [36J, 

Gm 2 Gm 2 \ / Gm\ Grrii 
ai • a 2 = I ai H — r — r • a 2 H — r 



Gmi Gm2 G 2 m\vri2 
- ai ■ — — r - a 2 ■ — ^-r H , 

which makes use of the lowest order equations of motion to get rid of the 3PN Lagrangian 
terms quadratic in the accelerations. This transformation should be iterated (as in the case 
of the a x • a 2 term) and combined with derivations by parts until only terms linear in the 
accelerations remain. 

Moreover, as terms like a x • a 2 appear also in the 2PN effective Lagrangian, the same 
technique must be applied at 2PN, but employing the equation of motion up to 1PN accuracy 
rather than just the Newtonian part written above; this gives a further contribution to the 
3PN Lagrangian, and all these terms nicely fit with all the previous ones to give eq. (174) 

of ID. 



IV. CONCLUSIONS 



We have computed the conservative dynamics of a gravitationally bound binary system 
within the framework of the PN approximation to general relativity to third post-Newtonian 
order. By a systematic use of the effective field theory methods for general relativity pro- 
posed by Goldberger and Rothstein it has been possible to automatize the computation 
and compute the effective Lagrangian by summing the contributions of eighty Feynman 
diagrams. We have computed the effective Lagrangian reproducing known results [3] and 
paving the way for the yet uncomputed 4PN dynamics. Beside the clear theoretical interest 
of this calculation, the 4PN Hamiltonian has a direct phenomenological impact as it enters 
the determination of the templates waveforms used in analyzing the data from the large 
interferometers LIGO and Virgo like the Effective One Body ones [28] . 
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